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Abstract 

Recently, it was observed that solutions of a large class of highly oscillatory second order linear ordi¬ 
nary differential equations can be approximated using nonoscillatory phase functions. In particular, 
under mild assumptions on the coefficients and wavenumber A of the equation, there exists a func¬ 
tion whose Fourier transform decays as exp(—//|£|) and which represents solutions of the differential 
equation with accuracy on the order of A” 1 exp(—pA). In this article, we establish an improved ex¬ 
istence theorem for nonoscillatory phase functions. Among other things, we show that solutions of 
second order linear ordinary differential equations can be represented with accuracy on the order of 
A -1 exp(— fi\) using functions in the space ,S’(M) of rapidly decaying Schwartz functions whose Fourier 
transforms are both exponentially decaying and compactly supported. These new observations play 
an important role in the analysis of a method for the numerical solution of second order ordinary 
differential equations whose running time is independent of the parameter A. This algorithm will be 
reported at a later date. 
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1. Introduction 


Given a differential equation 


y"(t) + \ 2 q(t)y(t) =0 for all a ^ t ^ b, 


( 1 ) 


where A is a real number and q : [0,1] —> R is smooth and strictly positive, a sufficiently smooth 
a : [a, b\ —> K is a phase function for (T]) if the pair of functions u, v defined by the formulas 

cos(a(t)) 


and 


u(t) = 


v(t) = 


k(i)| 1/2 

sin(a(t)) 


( 2 ) 


( 3 ) 


|a'(f )| 1/2 

form a basis in the space of solutions of £[]). Phase functions have been extensively studied: they were 
first introduced in |lj], play a key role in the theory of global transformations of ordinary differential 
equations , and are an important element in the theory of special functions [3,15, 6, 7|. 
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It was observed by E.E. Kummer in [1] that a is a phase function for 0 if and only if it satisfies the 
third order nonlinear differential equation 


(o/(t)) 2 = X 2 q(t) 


1 a'" {t) 

2 a'(t) 


4 \ a'{t) ) 


(4) 


on the interval [a, b]. The presence of quotients in © is often inconvenient, and we prefer the more 
tractable equation 


r "{t) - j ( r'(t )) 2 + 4A 2 (exp(r(f)) - q(t )) = 0 (5) 

obtained from ID) by letting 

a'{t) = A exp • ( 6 ) 

Of course, if r is a solution of 0 then the function a defined by the formula 

a{t) = A J exp du (7) 

is a solution of (0. We will refer to 0 as Kummer's equation and 0 as the logarithm form of Kum¬ 
mer 's equation. The form of these equations and the appearance of A in them suggests that their 
solutions will be oscillatory — and most of them are. However, there are several well-known exam¬ 
ples of second order ordinary differential equations which admit nonoscillatory phase functions. For 
example, the function 


a(t) = Aarccos(f), 


( 8 ) 


is a phase function for Chebyshev's equation 

y"{t) + + ^ !/(*) = 0 for a11 - 1 < t < 1. (9) 

Its existence is the basis of many numerical algorithms, including the fast Chebyshev transform (see, 
for instance, @])- Bessel's equation 

y"(t) + f 1 — ——\ y(t) = 0 for all 0 < t < oo (10) 


also admits a nonoscillatory p 
functions (see, for instance, |9|] 


hase function, although it cannot be expressed in terms of elementary 
)■ 


Exact solutions of 0 which are nonoscillatory need not exist in the general case. However, as we show 
in this article, when the coefficient q appearing in 0 is nonoscillatory, there exists a nonoscillatory 
function a such that 0, 0 approximate solutions of 0 in the space L 00 ([a, b ]) with accuracy on the 
order of A -1 exp(— y\). In order to make this statement rigorous, we will use the Fourier transform to 
quantify the notion of "nonoscillatory function." Accordingly, we assume that the coefficient q in 0 
extends to a strictly positive function on the entire real line. Moreover, we define the function x via 
the formula 

x{t) = f y/q{u) du , (11) 

Ja 

and let p(x) be twice the Schwarzian derivative of the variable t with respect to the variable x (see 
Section IZ8l) . We suppose that there exist constants T and // such that 

|p(£)l < r exp(-/i|£|) for all (eR, (12) 

that A > j l , and that A > 2T. Then there exists a function 5 in the space 5(E) of rapidly decaying 
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Schwartz functions (see Section l2TTb such that 


r{t) = log (q(t)) + 6(x(t)) (13) 

closely approximates a solution of ©, the support of the Fourier transform 5 of 5 is contained in the 
interval (—V2A, V2A), and 


£(£) < -r exp {—(jl |£|) for all |£| < V2A. 
A 


Moreover the function a defined via the formula 


a(t) = A f 

Ja 


exp 


r(u) 


du 


(14) 

(15) 


is a phase function for a second order differential equation of the form 

y"(t) + A 2 ^1 + j q{t)y{t) = 0 for all a ^ t 6, (16) 

where v is an element of S(M) whose (M) norm is on the order of exp (—//A). While a is not a phase 
function for the original equation ((Tj), the functions u, v obtained by inserting (flSl) into the formulas 
© and © approximate solutions of £Q with accuracy on the order of A 1 exp (—//A). 


The bound dl4l) implies that when S is approximated using various series expansions, the number of 
terms required to represent it to a specified precision is independent of A. For instance, the minimum 
number of terms in the Legendre expansion of the restriction of 5 to the interval [a, b] required to 
achieve a specified precision is independent of A. Assuming that q is nonoscillatory, it follows that r 
and a can be represented likewise; that is, using finite series expansions whose number of terms is in¬ 
dependent of A (it is in this sense that they are nonoscillatory). In other words: 0(1) terms are required 
to represent the function 6 which enables us to approximate solutions of CD) with accuracy on the order 
of O (A -1 ex p (— / /A)). This is in contrast to superasymptotic and hyperasymptotic expansions (see, for 
instance, [id. 11]) which approximate solutions of CD) to accuracy on the order of exp(— pA), but which 
require O(A) terms in order to do so. Note that we avoid discussing the Fourier transforms of r and a 
because they need not decay at infinity and so there is no assurance that their Fourier transforms are 
functions (as opposed to tempered distributions). 

The results presented in this article improve upon those of Hi which observed that solutions of (H]) 
can be represented with accuracy on the order of exp(—pA) using functions which are in L 2 (M) n 
C' 00 (M) and whose Fourier transforms decay exponentially. Here we show that the function 6 which 
represents the solutions of CD) is an element of the space S(R) of rapidly decaying Schwartz functions 
(see Section l2d]> and that its Fourier transform is both exponentially decaying and compactly sup¬ 
ported (so that 5 is entire). Moreover, we substantially reduce the constants appearing in the bounds 
on the decay of the Fourier transform of A and in the error of the associated approximations of the 
solutions of 0. Among other things, these new observations play an important role in the analysis 
of an algorithm for constructing a nonoscillatory solution a of 10) whose running time is independent 
of A. This algorithm allows for the numerical evaluation of solutions of second order linear ordinary 
differential equations of the form 0 using a number of operations which is independent of A. It will 
be reported at a later date. 


The remainder of this paper is organized as follows. Section [2] summarizes a number of well-known 
mathematical facts and establishes the notation which is used throughout this article. In Section0 we 
reformulate Kummer's equation as a nonlinear integral equation. The statement of the main result of 
this paper. Theorem fl2l is given in Section0and its proof is divided among Sections 0(6] and 0 
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2. Preliminaries 


2.1. Function spaces 


We denote by C (R) the set of continuous functions R —> C. If / e C (R) and, for each e > 0, there 
exists a compact set K such that |/(x)| < e for all x K, then we say that / vanishes at infinity. We 
denote the set of continuous functions which vanish at infinity by Cq (R). By C x (R) we mean the set 
of infinitely differentiable functions R —> C, and Cf (R) is the set of compact supported functions in 

C rao (R). 


We say that y e (7 00 (R) is a Schwartz function if y and all of its derivatives decay faster than any 
polynomial. That is, if 

sup < co (17) 

teK 


for all pairs i. j of nonnegative integers. The set of all Schwartz functions is denoted by S'(R); clearly, 
it contains the set Cf (R). We endow S’(R) with the topology generated by the family of seminorms 

k 

IMIfc = X! sup 

j =0 teR 




k = 0,1,2,... 


(18) 


so that a sequence {y n } of functions in S’(R) converges to y in S’(R) if and only if 


lim ||y n — y|L = 0 for all k = 0,1,2,... . 

n—>00 


(19) 


We denote the space of continuous linear functionals on S’(R), which are known as tempered distribu¬ 
tions, by S'^R). We endow S'{ R) with the weak-”' topology so that a sequence {co n } in S'( R) converges 
to co e 5 / (R) if and only if 

lim w n (y) = w(y) (20) 

n^co 

for all y e 5(R). We refer the reader to |iH] for a thorough discussion of the properties of Schwartz 
functions and tempered distributions. 


2.2. The Fourier transform 

We define the Fourier transform of a function / e S( R) via the formula 

^ rco 

7(0 = exp (-ixOf(x) dx. (21) 

J— 00 

The Fourier transform is an isomorphism 5(R) —> S(R) (meaning that it is a continuous, invertible 
mapping S( R) —> ,S’(R) whose inverse is also continuous). The formula 

(u),y) = (cj,y) (22) 

extends the Fourier transform to an isomorphism S"(R) —> S"(R). The definition (l22ll coincides with 
(I2B when / e L 1 (R). Moreover, when f e L 2 (R), 

C R 

/(£) = lim exp (—ix£)f(x) dx. (23) 

Owing to our choice of convention for the Fourier transform, 

T^g(o = hmo ( 24 ) 

and 

/ • g(£) = J ^ f Xf - v)div) dp (25) 
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whenever / and g are elements of L 1 (M). Moreover, 

1 


/(*) = 


2vr 


L 


exp(i®0/(0 


(26) 


whenever / and / are elements of L 1 (M). The observation that / is an entire function when / is a 
compactly supported distribution is one consequence of the well-known Paley-Wiener theorem. See 
I Idllla] for a thorough treatment of the Fourier transform. 


2.3. Convolution exponentials 
Formally, the Fourier transform of 


exp (/(s)) 


is the sum 


2 ^( 0 +m + ^ ^ * m 


(27) 

(28) 


2!(2vr) 3!(2vr) 2 

where ip denotes the Fourier transform of / and 6 is the delta distribution. The expression (l28l) . which 
is referred to as the convolution exponential of ip, is typically denoted by exp* [ip ]. 

We will not encounter the expression (l28l) in this article; however, we will consider the Fourier trans¬ 
forms of functions of the form 


exp(/(x)) - 1 


and 


exp(f(x)) - f(x) - 1. 

So, in analogy with the definition of exp*, we define exp* [ip] and exp| [V’j for ip e L 1 


formulas 


and 


expj [ip] (0 = ip(0 

ex P2 bP\ (0 = 


ip * ip(£) ip * ip * ip(£) 
2!(2vr) + 3!(2vr) 2 

ip * ip(fp) ip * ip * ip{£,) 


(29) 

(30) 
via the 

(31) 

(32) 


2!(2vr) 3!(2vr) 2 

That is, exp* [ip] is obtained by truncating the leading term of exp* [ip] and exp| [ip | is obtained by 
truncated the first two leading terms of exp* [ip]. By repeatedly applying the inequality 

11/ * slli < ll/IMMli) 

which can be found in lfl6h (for instance), we obtain 


Mi 


Ip * Ip 


ip * ip * ip 


00 

• < V 

Mi 

2!(2vr) 

1 

3!(2vr ) 2 

1 

- Zj 

n= 1 

n\(2'jr) n ~ 1 


(33) 

(34) 


and 


ip * ip 


ip * ip * ip 


. £ Mi . 

Mi / 

'Mi 

2!(2vr) 

1 

3!(2vr) 2 

1 

" — 2 n! ( 27r ) n_1 

4 7T 6XP l 

2vr 


(35) 


from which we see that the series (13Tb and d32|> converge absolutely in L 1 (M) when ip eM and therefore 
define L 1 (M) functions. Suppose that / is the inverse Fourier transform of ip e L 1 (M). For each 
nonnegative integer n, we define f n by 

{f(x)) k 


fn{x) = 


k =1 


k\ 


(36) 
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( 37 ) 


to exp(/(x)) — 1. Since 

ip * ip(0 ip * ip * 


HO 


We observe that {f n } converges in L 00 

fn{0 - HO . 2\{2ir) ' ' n!( 27 r) n_1 

for each nonnegative integer n, it follows from (l34l) that j/nj converges in if (R) to exp* [ip]- We 
conclude that the Fourier transform of exp(/(x)) — 1 is exp* [ip]- A nearly identical argument shows 
that the Fourier transform of exp(/(x)) — f(x) — 1 is exp| [ip]- 

Theorem 1. If f e S(M), then exp (f(x)) — 1 and exp (f(x)) — f(x) — 1 are elements of S(M). 


Proof. Since / e S(M), exp(|/(x)|) is bounded and 

k 


sup \x k f(x)\ < 00 
:ceR 


for any nonnegative integer k. Consequently, 

x k (exp (f(x)) - 1) 


sup 

xsR 


R sup 

x k f(x) 

seR 



(38) 


(39) 


for any nonnegative integer k. We conclude that exp (f(x)) — 1 decays faster than any polynomial. We 
observe that the n th derivative of exp(/(x)) — 1 is of the form 

P {}{x),f'[x),f"{x ),..., / (n) 0r)) exp {f{x)), (40) 

where P is a polynomial in n variables of the form 

P{xi,X 2 , ...,X n )= C k 1: k 2 ,-,k n x l lx 2 2 ■ ■ ■ X n n ■ ( 4i ) 

l^ki+k2 + ...+k rl ^n 

Since /, f, f", ..., /^ are elements of S(M), 

(/(x)) fcl (/'(x)) fc2 ...(/W(x)) fcri (42) 

decays faster than any polynomial whenever k \, & 2 , ■ ■ ■, k n are nonnegative integers not all of which 
are 0. We combine this observation with the fact that exp(/(x)) is bounded in order to conclude that 
the function ( f40l) decays faster than any polynomial; i.e., the n th derivative of exp(/(x)) — 1 decays 
faster than any polynomial. Therefore exp(/(x)) — 1 is in S'(M). Since exp(/(x)) — f(x) — 1 is obtained 
from exp(/(x)) — 1 by subtracting the Schwartz function f, it is also an element of ,S'(R). □ 

We combine Theorem |l] with the observation that the Fourier transform is a continuous mapping 
S(M) —» S(]SL) in order to obtain the following theorem. 

Theorem 2. If ip e S'(M), then exp* [ip\ and exp| [ip] are elements of S(M). 


2.4. The constant coefficient Helmholtz equation 

Under certain conditions on the function /, a solution of the inhomogeneous Flelmholtz equation 

y"{x ) + A 2 y(x) = f(x) for all x e M (43) 

can be obtained via the Fourier transform. For instance, the following theorem is a special case of a 
more general one which can be found in [17]. 


Theorem 3. Suppose that f e L 1 (M) n C(M), and that A is a positive real number. Then the function g defined 
by the formula 


1 f 00 

= 2\ J sm ( x \ x ~ y\) f(y) d y 


(44) 
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is twice continuously differentiable, 

g" (x) + X 2 g(x) 

mid 


g(0 


f(x) for all xel, 

m 

a 2 - e 2 ' 


(45) 

(46) 


We interpret the Fourier transform (l46l> of g as a tempered distribution defined via principal value 
integrals; that is to say that for all ip e S( R), 


m) 


.¥> ) = 



/(SM£) 


df — lim 
6—>0 


1 


/(£M0 

A + e 


dt, 


'h-A|>e A-£ e-ujj ?+A | >e 

We also note that the requirement that / is continuous ensures that (l45l) holds for all x e 
such an assumption on f, we are only guaranteed that (l45ll holds for almost all x e R. 


(47) 
i; without 


When / e L 2 (R) and A is real-valued, the integral (l44l defining the function g is not necessarily 
absolutely convergent and the expression 


m 

A 2 - £ 2 ’ 


(48) 


which is formally the Fourier transform of g, need not define a tempered distribution. If, however, the 
support of / is contained in (—A, A), then (l48l) is a compactly supported element of L 1 (R). In this case, 
we define g through the formula 


1 f 00 f(6) 

5(®) = ^:J exp(ix^) ^ df. (49) 

Since g is the inverse Fourier transform of a compactly supported function, it is entire. Moreover, the 
Fourier transform of 


g"(x) + A 2 g(x) 


(50) 


is 


from which we conclude that 


,2 _m_ , ,2 m 

? A 2 - £ 2 A 2 - £ 2 


m, 


(51) 


g"{x) + A 2 g{x) = f{x) (52) 

almost everywhere. Since both g and / are both entire, (l52l) in fact holds for all xel. We record these 
observations as follows. 


Theorem 4. Suppose that f e L 2 (R), that A is a positive real number, and that the support of f is contained 
in the interval (—A, A). Then the function g defined via the formula 

1 f 00 f(£) 

^ = 2tt J ex v( ix & x 2 - g 2 ^ ( 53 ) 

is entire, 

g"(x) + A 2 g(x) = f(x) for all x e R, (54) 

and 

m - iff?- ( 55 ) 


7 












The following variant of Theorem [3] can be found in 


Theorem 5. Suppose that f is continuous on the interval [a, b], and that A is a positive real number. Suppose 
also that y : [a, b] —> C is twice continuously differentiable, and that 

y"(x) + A 2 y(x) = f(x) for all a ^ x < b. (56) 


Then 


y{x) = y{a ) + y'(a)(x 



sin (A (x 


u)) f(u) du for all a ^ x < b. 


(57) 


2.5. Modified Bessel functions 


The modified Bessel function K v (t) of the first kind of order v is defined for t e M and u e Cby the 
formula 


r 00 

K v (t) = exp (—t cosh (t)) cosh(z^f) dt. 

Jo 

The following bound on the ratio of K v+ i to I\ u can be found in Ml 


(58) 


Theorem 6. Suppose that t > 0 and v > 0 are real numbers. Then 

K u+ i(t) v + y/u 2 + t 2 2 v 

K v (t) < t " t + 


(59) 


2.6. The binomial theorem 

A proof of the following can be found in l(2otl . as well as many other sources. 


Theorem 7. Suppose that r is a real number, and that y is a real number such that \y\ < 1. Then 


(i + y) r 


y F(r + 1) fc 

f^nk + mr-k + iy ■ 


(60) 


2.7. Frechet derivatives and the contraction mapping principle 


Given Banach spaces X, Y and a mapping / : X —> Y between them, we say that / is Frechet differ 
entiable at x e X if there exists a bounded linear operator X —> Y, denoted by f x , such that 

II f(x + h) - f(x) - /' [h] 


lim ■ 

h—>0 


\h\\ 


= 0 . 


(61) 


Theorem 8. Suppose that X and Y are a Banach spaces and that f : X —► Y is Frechet differentiable at every 
point of X. Suppose also that D is a convex subset of X, and that there exists a real number M > 0 such that 

ll/ill < M (62) 


for all x £ D. Then 


f(x) ~ f(y )|| M\\x - y 


(63) 


for all x and y in D. 
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Suppose that / : X —*• X is a mapping of the Banach space X into itself. We say that / is contractive 
on a subset D of X if there exists a real number 0 < a < 1 such that 

11/0*0 -/(y)||< &\\x-y\\ (64) 

for all x,y e D. Moreover, we say that {x n }® =0 is a sequence of fixed point iterates for / if x n+ \ = f(x n ) 
for all n ^ 0. 

Theorem [8] is often used to show that a mapping is contractive so that the following result can be 
applied. 


Theorem 9. (The Contraction Mapping Principle) Suppose that D is a closed subset of a Banach space X. 
Suppose also that f : X —*• X is contractive on D and f(D ) c: D. Then the equation 

x = f(x) (65) 


has a unique solution a* e D. Moreover, any sequence of fixed point iterates for the function f which contains 
an element in D converges to a*. 


A discussion of Frechet derivatives and proofs of Theorems [8] and [9] can be found in, for 


instance, [21], 


2.8. Schwarzian derivatives 


The Schwarzian derivative of a smooth function / : R —> R is 

nt) 3 (f"(t )\ 2 

17,1 f'{t) 2 \f'(t)) 


( 66 ) 


If the function x(t ) is a diffeomorphism of the real line (that is, a smooth, invertible mapping M —> M), 
then the Schwarzian derivative of x(t) can is related to the Schwarzian derivative of its inverse t(x); 
in particular, 

{x,t} = ~ ( 67 ) 

The identity (1671) can be found, for instance, in Section 1.13 of |[f«]. 


2.9. A bump function 


It is well known that the function defined by the formula 

- 1 K 


" (L exp (^3i) du ) £„ exp (^i) 


( 68 ) 


is an element of C 00 (M) such that <p(£) = 0 for all ^ < — 1, p(^) = 1 for all £ > 1 and 0 ^ <p(£) ^ 1 
for all ^ e R (see, for instance, fhil). We suppose that A is a positive real number and define the bump 
function b via the formula 


where 


and 


m 



a = 



V2X + A 


c — A \/2A — A 
4 = 4 



(69) 

(70) 

(71) 
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We observe that b is an element of C® (M), the support of b is contained in (- v'2A, 0 ^ b(f) ^ 1 

for all (eK, and b(f) = 1 for all / ^ A. Since b is an element of (7® (M), its inverse Fourier transform 
b is an entire function and an element of 


2.10. The Liouville-Green transform 


The Liouville-Green transform is a well-known tool for analyzing the variable coefficient Helmholtz 
equation 


y"(t) + A 2 q{t)y{t) = /(f). 


The following can be found in Chapter 2 of [22] (for instance). 


(72) 


Theorem 10. Suppose that q : [a, 6] —> M is twice continuously differentiable and strictly positive, that 
f : [a, 6] —> C is continuous, that the function x is defined by the formula 

x(t) = f V q(u ) du, (73) 

Ja 

and that the function p is defined by the formula 


p{t) 


_J_ (5 (q'jt)\ 2 q"jt)\ 

q(t) \ 4 V ?(*)) q(t) ) ' 


(74) 


Suppose also that y : [a, 6] —> R is twice continuously differentiable, and that 

y''{t) + \ 2 q{t)y{t) = /(f) for all a ^ f ^b. 


(75) 


Then the inverse t(x) of the function x{t) is continuously differentiable, and the function <p : [0, m(6)] —> M 
defined by the formula 

(p{x) = (g(f(x))) 1/4 y(f(x)) (76) 


is the unique solution of the initial value problem 

ip"(x) + \ 2 ip(x) = (g(x))~ 3 /, 4 /(x) — ^p(x)<p(x) for all 0 ^ x ^ x(b) 
' <p(0) = (<?(a)) 1/4 ?/(a) 

<p'(0) = q'{a)(q{a)y 5/4 y{a) + (<?(a)) _1/ V (a). 


Remark 1. We observe that due to Formula (1671 ). the function p(x) appearing in ([771 is twice the Schwarzian 
derivative of the inverse f(x) of the function x(t) defined in (f73l >. That is, p(x) = 2{t, x}. 


2.11. Gronwall's inequality 

The following well-known inequality can be found in, for instance, [23]. 

Theorem 11. Suppose that f and g are continuous functions on the interval [a, b] such that 


/(f) ^ 0 and gift ) ^ 0 for all a ^ f ^ b. 

Suppose further that there exists a real number C > 0 such that 

(78) 

/(f) ^ C + f f{s)g(s) ds for all a ^ f ^b. 

J cl 

(79) 

Then 


/(f) ^ C exp ^J g{s) ds^j for all a < f < b. 

(80) 
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3. Integral equation formulation 


In this section, we reformulate Kummer's equation 


(«'(*)) 2 = A 2 g(t) 


1 a"'(t) 

2 a/(f) 


3 

4 V «'(*)/ 


( 81 ) 


as a nonlinear integral equation. We assume that the function q has been extended to the real line and 
we seek a function a which satisfies dHTb on the real line. 


By letting 

[a(t)Y = A 2 exp (r(t)) (82) 

in (f8Tb . we obtain the equation 

r"(t) — i (/(f)) 2 + 4A 2 (exp(r(f)) — q(t)) = 0 for all t e R. (83) 

We next take r to be of the form 


r(t) = log(g(f)) +S(t), 


which results in 


1 v'(t) , 


$"{t) - ~ A (exp(5(f)) - 1) = q(t)p(t), for all t 


2 q(t) ' " 4 

where p is defined by the formula 


Pif) = 


J_( 5 /VEA 2 q"jt) 

q(t) \ 4 V q(t)) q(t) 


Expanding the exponential in a power series and rearranging terms yields the equation 

«"(<) - + 4A 2 «(t)<5(*) - 1 (4'{*)) 2 + 4A 2 ,(t) ( Ml! + Ml! + ...)= 9 (t)j>(t) 


2 q(t) 

Applying the change of variables 

transforms d87l) into 


*(*) = f Voiuj 

Ja 


u ) du 


S"(x) + 4A 2 6 (x) - - (<5'(x)) 2 + 4A 2 


(5(x)) 2 , (,5(x)) ; 


3! 


= p(x) for all lel 


(84) 

(85) 

( 86 ) 

(87) 

( 88 ) 

(89) 


At first glance, the relationship between the function p(x) appearing in (l89l) and the coefficient q(t) in 
the ordinary differential equation ([TJ) is complex. However, the function p(t) defined via (l86l) is related 
to the Schwarzian derivative (see Section fZ8l) of the function x(t) defined in (l88l) via the formula 

pit) = ~q(t) { x,t i = ~ 2 (ytx) ( 9 °) 

It follows from l90l) and Formula (l67l) in Section f2.81 that 

p(x) = 2 {t, x}. (91) 

That is to say: p, when viewed as a function of x, is simply twice the Schwarzian derivative of t with 
respect to x. 
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It is also notable that the part of (I89l> which is linear in 5 is the constant coefficient Helmholtz equation. 
This suggests that we form an integral equation for (f89l) using a Green's function for the Helmholtz 
equation. To that end, we define the linear integral operator T for functions /ef 1 (R) n C (R) via the 
formula 

T [/] (x) = ^ J sin (2A \x — y\) f(y) dy. (92) 

whose Fourier transforms have support in 


We extend the domain of T to include functions / e L 2 
the interval (—2A, 2A) through the formula 


T 




4A 2 — £ 2 

Introducing the representation 

5{x) = T[a\ (x) 

into 1891) yields the nonlinear integral equation 

a{x) = S [T [cr]] (x) + p(x) for all x e R, 
where S is the nonlinear differential operator defined by the formula 

(/M ) 2 . (/M ) 3 , (/W)‘ 


s[/]W-^^ML-4A 2 


2 ! 


3! 


4! 


(93) 

(94) 

(95) 

(96) 


According to Theorems |3] and |H if a is a solution of the integral equation (l95l) and either a e L 1 (R) n 
C (R) or a e L 2 (R) and the support of a is contained in (—2A, 2A), then the function 6 defined via 
formula l94l) is a solution of <f89l) . Moreover, the function r defined via the formula 


is a solution of (1831) . and 


is a phase function for dTJ. 


r(t) = log(g(t)) + «5(x(t)) 
' r(u) 


a(t ) = A exp 

Ja 


du 


(97) 

(98) 


4. Existence of nonoscillatory phase functions 


The nonlinear integral equation (l95l> is not solvable for arbitrary p. However, when the Fourier trans¬ 
form of the function p decays exponentially, there exists a function a whose Fourier transform is com¬ 
pactly supported and a function v of magnitude on the order of exp(—pA), where p is a real constant, 
such that 

a{x) = S [T [cr] ] (x) + p(x) + v(x) for all x 6 R. (99) 

The following theorem, which is the principal result of this article, makes these statements precise. Its 
proof is given in Sections |5j [6] and [7| 


Theorem 12. Suppose that q e C 00 (R) is strictly positive , that x(t) is defined by the formula 

*(*) = f Voi 1 

Jo 

and that the function p defined via the formula 


(u) du , 


p{x) = 2 {t,x} 


( 100 ) 

( 101 ) 
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is an element ofS(M). Suppose furthermore that there exist positive real numbers X, T and p such that 


and 


A > 2 max ( —, T 
l 1 


( 102 ) 


\p{£,)\ ^ r exp (—p |£|) for all (eM. (103) 

Then there exist functions a and v in S(R) such that a is a solution of the nonlinear integral equation 

(t(x) = S [T [cr]] (x) + p(x) + u(x), for all x e R, (104) 

1^(01 < + 


2T 


Texp(— p\£\) for all |£| ^ V2X, 


(105) 


and 


d(f) = 0 for all |£| > \/2A, 


v 




exp (—pX). 


Remark 2. By combining (I102D zvith (I105P . we see that 

|ct(£)| sS 2r exp (-/u|f |) 
for all |£| ^ V2A. Similarly, we conclude from (11021 ) and (12071 ) that 

3F 

Halloo < — exp(-pA) • 


(106) 

(107) 


(108) 

(109) 


Suppose that a and v are the functions obtained by invoking Theorem [T2l and that x(t) is the function 


defined by the formula 




x(t) = V Q( u ) du. 

(110) 

We define 5 by the formula 




5{x) = T [cr] (x), 

(111) 

r by the formula 




r(t) = log (q(t)) + 6(x(t)), 

(112) 

and a by the formula 




du. 

(113) 


From the discussion in Section |3j we conclude that d(x) is a solution of the nonlinear differential 
equation 

8"(x) + 4A 2 5(x) = S [5] (x) + p(x) + u(x) for all ieI, (114) 


that r(t) is a solution of the nonlinear differential equation 


r"(t ) — ~(r'(t)) 2 + 4A 2 (exp (r(t)) — q(t )) = q(t)v(t) for all t e 


and that a is a solution of the nonlinear differential equation 


(aft)) 2 = A 2 


v(tf 

4A 2 


, , 1 a’" {f) 

**> - iiM 


q"(Q 

a'{t) 


for all t e 


(115) 


(116) 
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From (11161) . we see that a is a phase function for the second order linear ordinary differential equation 

y"(t) + A 2 ^1 + q(t)y(t) = 0 for all a ^ f ^ b. (117) 

Since the magnitude of v is small, we expect that the difference between solutions of (11171) and those of 
dTj) will be small as well. Indeed, the following theorem follows easily by applying the Liouville-Green 
transform and invoking Gronwall's inequality. 


Theorem 13. Suppose that q is continuous and strictly positive on the interval [a, b\, that f is a continu¬ 
ous function [a,b] —> C, and that A is a positive real number. Suppose also that z : [a, 6] —*• C is a twice 
continuously differentiable function such that 

z"(t) + X 2 q(t)z{t) = 0 for all a ^ t ^ b, (118) 

and that y : [a, b] —» C is a twice continuously differentiable function such that 

y"(t ) + X 2 q(t)y(t) = /(f) for all a t ^ b, (119) 

y(a) = z(a), and y'{a) = z'{a). Then there exists a positive real number C such that 

\y(t) — z(t)\ ^ y sup |/(t)| for all a ^ t ^ b. (120) 

The constant C depends on q but not on X or f. 


Proof. We define the function A by the formula 

A(t) = y{t) - z(t) 

and observe that A is the unique solution of the boundary value problem 

( A "(t) + X 2 q(t)A(t) = /(t) for all a ^ f ^ b 
"I A(a) = A/a) = 0. 

We define the function x(t ) via the formula 

x(t) = V q(u) du 

Ja 

and use t(x) to denote its inverse; moreover, we define the function p(t) via the formula 


p(t) = 

PU q(t)U{q(t)) q(t ) 


q(t) V q{t) J q(t ) J ' 

According to Theorem [lOj the function defined via 

¥>(*) = (^(f(^))) 1/4 A(i(®)), 

is the unique solution of the initial value problem 

<p"(x) + X 2 p(x) = ( q(x )) 3 ^ 4 f(x) — -p(x)<p(x) for all 0 ^ t ^ x(b) 

<p( 0) = <p'(0) = 0. 

We conclude from (11261) and Theorem [5] that 

= j J| sin (A(x-u)) ((q(u))~ 3 / 4 f(u) - ^p(u)<p(u)^j du 

for all 0 ^ x ^ x(b). We let 

Ci = x[b) ■ sup (q(u))~ 3 ^ 
a^u^b 


( 121 ) 

( 122 ) 

(123) 

(124) 

(125) 

(126) 

(127) 

(128) 
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and 


11 oo = sup 1 / 0)1 
a^u^b 


We observe that 


— f sin (A(x — u)) ( q(u )) 3 ^ 4 f(u ) du 

* Jo 


^ ^ ll/lloo 


for all 0 x ^ x(b). Now we let 


C 2 = sup p(u) 

a^u^b 


and observe that 


i r @2 r 

4aJ 0 sin ( A ( x ~ u )) P( u )<t>( u ) du ^ 00)1 du 


for all 0 < x < x(b). We combine (11271) . (11301) and (11321) in order to conclude that 

00)1 < ^-ll/lloo + 77 J 00)1 du 


(129) 

(130) 

(131) 

(132) 

(133) 


A 4A 

for all 0 ^ x ^ x(b). By invoking Gronwall's inequality (which is Theorem [Tl] in Section 12.1 1 1 ) we 
conclude that 


00)1 < — ll/lloo exp 


C 2 x 

4A 


for all 0 ^ x ^ x(b). We conclude from (11211) . (1125b and (11341) that 

I y(t) - z(t)\ ^ jll/Hoo 

for all a ^ t ^ b, where C is defined via the formula 

C = C X exp ( C2 ^ b) ) ■ sup {q(u))~ 1/A . 
V / a^u^b 


(134) 

(135) 

(136) 

□ 


By applying Theorem[l3]to dll7l) and ([TJ we obtain the following. 

Theorem 14. Suppose that the hypotheses of Theorem H2l are satisfied, that a and v are the functions obtained 
by invoking it. Suppose also that a is defined as in ill3t . and that u, v are the functions defined via the formulas 

cos 00) 


and 


u{t) = 

v{t) = 


V a '(*) 

sin(a(t)) 


Then there exist a constant C and a basis {«, u} in the space of solutions of {IJ such that 


and 


(j 

\u(t) — u(t)\ ^ — exp (— pX) for all a ^ t ^ b 
A 


c 

\v(t) — v(t)\ ^ — exp (— pX) for all a^t^b. 
A 


(137) 

(138) 

(139) 

(140) 


The constant C depends on the coefficient q appearing in ([/]), but not on the parameter A. 

The rest of this article is devoted to the proof of Theorem [121 It is divided among Sections OHI and 0 
The principal difficulty lies in constructing a function v such that (l99l) admits a solution. We accom- 
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( 141 ) 


plish this by introducing a modified integral equation 

(T b {x) = S[T b [cr 6 ] ] (x) + p(x), 

where T b is a "band-limited" version of T. That is, T b [f] is defined via the formula 

n\7m = T\fKmo, M 

where &(£) is the Gf (E) bump function given by Formula d69l ) of Section [2791 In Section [21 we apply 
the Fourier transform to dl41b and use the contraction mapping principle to show that under mild 
conditions on p and A the resulting equation admits a solution. This gives rise to a solution a b of dl41ll . 
In Section [6} we show that if the Fourier transform of the function p is exponentially decaying, then 
the Fourier transform of cr b is as well. In Section[7]we use the solution a b of dl41l) in order to construct 
functions a and v which satisfy d99l) . Moreover, we show that a can be taken to be an element of the 
space 5(E) of rapidly decaying Schwartz functions (see Section l2Tl) . that the Fourier transform of a 
is compactly supported and exponentially decaying, and that v is an element of 5(E) whose L 00 (E) 
norm decays exponentially with A. 


5. Band-limited integral equation 


In this section, we introduce a "band-limited" version of the operator T , use it to form an alternative 
to the integral equation (l95l) . and apply the contraction mapping principle in order to show that this 
alternate equation admits a solution under mild conditions on the function p and the parameter A. 

Let b(fi) be the bump function defined via formula (l69l) so that 

1. S(£) = 1 for all |£| < A, 

2. 0 < b(fi) < 1 for all (el, and 

3. b is supported on a proper subset of (— y/2X, \/2A). 

We define the operator T b [f] for functions / e Co (E) such that / el 1 (E) via the formula 

= (m 

We will refer to T b as the band-limited version of the operator T and and we call the nonlinear integral 
equation 

(j b (x) = S [T b [<7b]] (x) + p(x) for all xel (144) 

obtained by replacing T with T b in d95l) the "band-limited" version of (|95ll . 


Remark 3. The function 


m 

4A 2 - C 2 


(145) 


is an element of Cf (E) since the support ofb is bounded away from the points +2A at which the denominator 
vanishes. Consequently, (11431 ) is a compactly supported tempered distribution and T b [/] is an entire function 
whenever f is a tempered distribution. For our purposes, it suffices to know that T b [f] is defined for f e L 1 (E). 


It is convenient to analyze (|144l I in the Fourier domain rather than the space domain. We denote by 
W b and XV b the linear operators defined for f e if (E) via the formulas 


w b [f](0 


f(0 


HO 

4A 2 - C 2 


(146) 
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and 


w b [fm = m-£0^, (M7) 

where b(£) is the function used to define the operator T b . We define functions ip(£) and w(£) using the 
formulas 

= <?&(£) (148) 

and 

w{£,)=P(0- (149) 

Finally, we define R [f] for functions / e Jj (R) via 

R [/] (£) = [f] * W b [f] (0 - 4A 2 expl [W b [f] ] (0 + w(Z), (150) 

where exp| is the operator defined by Formula (f32l) of Section [231 Applying Fourier transform to both 
sides of (11441) results in the nonlinear equation 

m = R&no- (i5i) 

The following theorem gives conditions under which the sequence {ip n }fp = 0 of fixed point iterates for 
(11511) obtained by using the function w defined by (11491) as an initial approximation converges. More 
explicitly, ip o is defined by the formula 

i’oiO = w(0> (452) 

and for each integer n ^ 0, ip n + 1 is obtained from ip n via 

lP n+ l(0=R[^n](0- (153) 

Theorem 15. Suppose that A > 0 is a real number, and that w is an element of L 1 (R) such that 

IM|i ^A 2 . (154) 

Then the sequence {ip n } defined by (I152D and (I153D converges in L 1 (R) norm to a function ip e L 1 (R) such 
that 

ip(£) = R [ip] (0 for all (eK. (155) 


Proof. We observe that the Frechet derivative (see Section 12.71) of R at / is the linear operator R'j : 
L 1 (R) —> L l (R) given by the formula 


R'f [h] (0 = 


W b [f]*w b [h] (0 

47r 


- 4A 2 exp* 


W b [f] 

2vr 


*W b [h]{i). 


From formulas (11461) and (11471) and the definition of &(£) we see that 

ll/ll 


and 


W b [f] ^ 11/1,1 


i V2A 

for all / e L 1 (R). We combine dl56b with (l34b in order to conclude that 

||W 6 [/]||i 


1 2A 2 

\\ R 'f Will ^ w b[f] W b [h] + - \W b [/] ||r exp 

J 1 47T 1 17r 


2t r 


|W 6 [/i]||i 


(156) 

(157) 

(158) 

(159) 
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for all / and h in L 1 (R). By inserting (11571) and dl58b into (11591) we see that 


KMIIi 


< 


H/llillhjli , 2A 2 ||/|ji \\h\\ 


< 


8ir\ 2 

/Mi 

\8n\ 2 


7r 2A 2 2A 2 

\f\\i 


exp 


(Ml d 

\4:TTX 2 J 


2nX 2 


exp 


/Mi 

V4vrA 2 


( 160 ) 


for all / and h in L 1 


. Similarly, by combining (11501) . (l35ll . (1157b and (1158b we conclude that 

W b [f] Hi' 


ll *[/]| li ^|^[/]'' 2 


— II W b [/]||i exp 
7r 


2t r 


+ \\ w h 




ll/Ilf 


I/ll? 


exp 


(Mk 

\.47rA 2 


(161) 


+ \\w 1 


167tA 2 47tA 2 

whenever / e L 1 (M). We now let r = nX 2 and denote by ft the closed ball of radius r centered at 0 in 
L 1 (R). Suppose that / £ L 1 (R) such that 

\\f\\l < r = ttA 2 , 

and that 

"" 7tA 2 


II' 1 ’111 < 2 


(162) 

(163) 


We insert (11621 ) and (11631 ) into dl61b in order to obtain 


IW]l!i 




167tA 2 
r 


47tA 2 

r 


exp 


167tA 2 47tA 2 


exp 


(Uh\ + 

\4ttA 2 / 

(in l\ 

\47rA 2 J 


— r 


1 


1 


= 'l6 + 4 eXP '- 

9 

y — r, 

10 


\ I r 


(164) 


from which we conclude that R maps Q into itself. Next, we insert (11621) into (11601) in order to obtain 


II*/ [h]l 




< 


8ttX 2 

1 1 


27tA 2 


exp 


■ exp 


47tA 2 

h\\i 


\h\\ 


(165) 


< — \\h\\i, 

10 11 11 


which shows that R is a contraction on Q. We now invoke the contraction mapping theorem (The¬ 
orem [9] in Section 12.71) in order to conclude that any sequence of fixed point iterates for (11511) which 
originates in 0 will converge in L 1 (R) to a solution of (11511) . Since {d ? ,} is such a sequence, it converges 
in L 1 (R) to a function ip such that 

= R \ip] (£) for almost all (eM. (166) 

We change the values of ip on a set of measure zero (which does not affect R [ip]) in order to ensure 
that (11661 1 holds for all (eR. □ 


6. Fourier estimate 

In this section, we derive a pointwise estimate on the solution ip of Equation d!51b under additional 
assumptions on the function w. 
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Lemma 1. Suppose that p and C are real numbers such that 

0 C < p. 

Suppose also that f e L 1 (M), and that 


Then 


lexps [/] (01 < 7r~ exp(-H^I) 1 + ^ exp 
Z7T jJj 


c 

2irp 

where exp| is the operator defined in d32l >. 

Proof. We let 

3i (0 = Cexp(-/x|£|) 

and for each integer m > 0, we define g m +1 in terms of g m through the formula 

3m+l(0 = 7^ 9m * 3l(0- 

We observe that for each integer m > 0 and all £ e E, 

'C\f \\ m ~ 1/2 K m - 1 / 2 (p\f\) 



(167) 

(el 

(168) 

ex pf^l£l) f oral1 £ e]R ’ 

(169) 


3m (0 = 2 VhC 


2tt 


T(m) 


(170) 

(171) 

(172) 


where K v denotes the modified Bessel function of the second kind of order u (see Section [231) . By 
repeatedly applying Theorem[6]of Section [2.51 we conclude that for all integers m > 0 and all real t, 

K m . 1/2 (t)«K V2 ( t ) n 


- *■/><*> (?) 


m_1 r (±±* + m _ i) 

rT¥) 


(173) 


We insert the identity 

k i/ 2 (t) = A^exp(-t) 

into dl73l) in order to conclude that for all integers m > 0 and all real numbers t > 0, 


K m — 1/2 W 


< 


0F / A 1/2-m r(±±* + m-l) 


2 V 2 


exp(— t) 


By combining dl75l) and dl72l) we conclude that 


3m(0 E C exp(—p|£|) ( — 


i-i r 


r(¥) 

i+mICI 


+ m — 1 


r(m)r 


i+Mki 

2 


(174) 


(175) 


(176) 


for all integers m > 0 and all £ E 0. Moreover, the limit as £ —» 0 of each side of (11761) is finite and the 
two limits are equal, so (11761) in fact holds for all £ e R. We sum (11761) over m = 2,3,... in order to 
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conclude that 


expa [ 9 ] (0 < C'exp(-^l^l) 

m=2 '^9/ 


m— 1 | 




+ m — 1 


T{m + l)r(m)r 


c 


1 +A t l'g| 


(177) 


Cexp(-/x|£|) — . / 

,^1 W/ r(m + 2)r(m + i)r( 


1+/4SI 

2 


for all (eR. Now we observe that 

1 


T(m + 2) 

Inserting dl78l) into (1177b yields 




for all m = 0 , 1 , 2 ,.... 


expl 


[ 9 ] (0 < C'exp(-^l) 2 


C 


Y I 1 +A t l?l 


m 


(178) 


(179) 


Tr'i V 2 W r(m + i)r(^) 

for all (eR. Now we apply the binomial theorem (Theorem[7]of Section [2.61 ). which is justified since 
C < a < 2irp, to conclude that 


l + /x|g| 

expl [g] (0 < Cexp(-p|C|) 1(1- I - 1 


= Cexp(—p|£|) 


^exp ^ 


2-JTfl 

1 + 


(180) 


log 


1 - 


- 1 


for all (et. We observe that 


exp(.x) — 1 ^ x exp(x) for all x ^ 0 , 


and 


0 ^ log 


1 — x 


^ 2x for all 0 ^ x ^ 


2vr 


(181) 

(182) 


By combining dl81l) and d!82l) with d!80b we conclude that 


ex P 2 [9] (0 < Cexp(-/x|f|) 
C 2 

< ^: ex P(-^l?l) 


l + hl£l 

2 

1 + mI£I 

9 


log 


1 - 


c 

27 Tfl t 


exp 


1 + a 4 !^! 


log 


1 - 


c 

2l TfL t 


exp 


c 

27 T/JL 


exp ( 9-\(\ 


(183) 


for all (el. Note that in dl 83b , we used the assumption that C < a in order to apply the inequality 
d!82l) . Owing to d! 68 l) . 

l ex P 2 [/] (01 < ex P 2 [9] (0 for all (e1 (184) 

By combining this observation with d!83l) . we obtain d!69b . which completes the proof. □ 

Remark 4. Kummer's confluent hypergeometric function M(a, b, z) is defined by the series 


M(a, M ) , 1 + ^ + <“ )2 " 2 ■ {ahz3 


b ( 6)2 2! (6) 3 3! 


where (a) n is the Pochhammer symbol 


/ \ r(a + n) , , , . 

(«)n = —prry— = a(a + l)(a + 2 )... (a + n - 1 ). 
T(a) 


(185) 


(186) 
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By comparing the definition of M(a, b, z) with (1 1771 ). we conclude that 

l ex P2 [/] (61 < C'exp(- Ai |e|) (m ^11^1,2, - l) for all £ e R (187) 

provided 


|/(6I < C'exp(-//|^|) for all (eR. (188) 

The weaker bound (I169D is sufficient for our immediate purposes, but formula (11871 ) might serve as a basis for 
improved estimates on solutions of Kummer's equation. 


The following lemma is a special case of Formula (11721) . 


Lemma 2. Suppose that C ^ 0 and p > 0 are real numbers, and that / el 1 (R) such that 


|/(6I ^ C exp (~p\£\) for all (eK. 

Then 

|/*/(6l<C' 2 exp(-H6)(^f^) forall (eR. 


(189) 

(190) 


We will also make use of the following elementary observation. 

Lemma 3. Suppose that p > 0 is a real number. Then 

exp(-/x|6)|6 forall (eR. (191) 

p exp(l) 


We combine Lemmas Q] and |2] with dl57l) and d l 581 ) in order to obtain the following key estimate. 


Theorem 16. Suppose that T>0,A>0, p> 0 and C > 0 are real numbers such that 

0 ^ C < 2pX 2 . 

Suppose also that f e L 1 (R) such that 

1/(61 < Cexp(-p\£\) forall |£| sc V2A, 

and that w e L 1 (R) such that 

M6I ^ Texp(— p\£\) forall (eR. 

Suppose further that R is the operator defined via dl50l >. Then 


\R[f] (61 < exp(-p|^|) 


(J_ 1 


A 2 


p J ^16^ ' 27t ^ \AttX 2 pJ 1 ^4^ A 2 


C 


exp 


(■ 


C 


■16 


for all (el. 


(192) 

(193) 

(194) 

(195) 


Proof. We define the operator R\ via the formula 


Ri if] (0 = ^w b [f]*w b [f](0 

(196) 

and Ro by the formula 


R 2 [f] (6 = —4A 2 ex P 2 \W b [/]] (6> 

(197) 

where W b and W b are defined as in Section [5] Then 


R[f] (6 = Ri [/] (6 + ^2 [/] (6 + ^(6 

(198) 
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for all £ e R. We observe that 


Wbimt) 


< 


c 

V2\ 


exp(—p|£|) for all £ e 


By combining Lemma [2] with (119911 we obtain 


\Ri [/] (0\< 


c 2 


167tA 2 


exp(-q|£|) 


1 + h-\£\ 


for all £ e 


Now we observe that 


(199) 


( 200 ) 


\W b [f] (£)| ^ ^exp(-/i|£|) for all £ e R. (201) 

Combining Lemma [T| with (12011) yields 

\ R i [/] (01 « OXP (i£^) exp (j^lfl) (202) 

for all £ e R. Note that (11921) ensures that the hypothesis ( 11671) in Lemma [[] is satisfied. We combine 
(12001) with d202l) and dl94l) in order to obtain d!95l) . and by so doing we complete the proof. □ 

Remark 5. Note that TheoremWfUonly requires that /(£) satisfy a bound on the interval [—a/2A, \/2A] and not 
on the entire real line. 


In the following theorem, we use Theorem H6l to bound the solution of (11511) under an assumption on 
the decay of w. 


Theorem 17. Suppose that A > 0, p > 0 and T > 0 are real numbers such that 

1 


A > 2 max < T 




Suppose also that w e L 1 (R) such that 

|w(£)| ^ Texp(—/r|£|) for all £ e R. 

Then there exists a solution £) of U51t in L 1 (R) such that 

|^(£)| < ^1 + ^0 Texp(—/x|£|) for all |£| ^ V2\ 


and 


1^(01 < (l + f) re xp(-(^-^x)l£l) M all £e 


(203) 

(204) 

(205) 

(206) 


Proof. From d203b and d204b we obtain 

fo° T A 2 

Mil < r exp(-// |£|) d£ = — < —. (207) 

J-00 h 4 

It follows from Theorem [L5l and d207D that a solution £>(£) of d!51b is obtained as the limit of the se¬ 
quence of fixed point iterates Mn(£)} defined by the formula 

MO = (208) 

and the recurrence 

^ n+ i(£) = i2[^ n ](£). (209) 

We now derive pointwise estimates on the iterates ip n {0 in order to establish d205b and d206l) . 
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We denote by {/3k} be the sequence of real numbers generated by the recurrence relation 

Am = |j + r (210) 

with the initial value 


Po = r. (211) 

From mathematical induction and (12031) . we conclude that {Pk } is a monotonically increasing sequence 
which converges to 

13= X- VA 2 - 2rA. (212) 


We also observe that 


Pn < P < 



r < 2 r. 


(213) 


Now suppose that n ^ 0 is an integer, and that 

IVUOI < Ai exp(-/x|f |) for all |£| ^ V2A. (214) 

When n = 0, this is simply the assumption (12041) . The function ^n+i(?) is obtained from d n (0 via the 
formula 


Vwi (Z) = RW](0- 

We combine Theorem [16] with (1215b and (12141) to conclude that 


|Y’n+l(£)l < exp(-//|£|) 


A 2 


( 




( exD (A 


U^ + ^ exp v«v exp u^ 


M\ 


for all £ £ M. The hypothesis dl92l) of Theorem Il 6 l is satisfied since 

p n < 2T ^ 2A 2 n 

for all integers n ^ 0. We restrict £ to the interval [—-\/2A, V2A] in (1216b and use the fact that 

1 1 

/xA < 2’ 

which is a consequence of (12031) . in order to conclude that 


|^n+l(OI < exp(-//|£|) 
< exp(-//|£|) 


Pi (l + vA\\ ( 1 , 1 ( /?„ \ ( Pn 

~ 1 ^^) \ 167T + 2ir 6XP ) exp \47 tA 2 


A 2 

A 


Vzx 




r * 


i^ + ^ exp ^ exp (ivsi; 


f ^ \ 


for all |^| ^ \/2A. By combining d219l) with inequality 


Pn 2T 

— ^ ^ 1 for all n ^ 0 

A A 


and the observation that 




(l 1 / 1 

I- 1 -exp l — 


\167r 27 t 


we arrive at the inequality 




8n 


exp 


2V2vr 


1 

< 2 ’ 


\^n+l(0\ ^ ( Ta + r ) 


= Pn +1 exp( —/i|^|), 

which holds for all |£| ^ y/2X. We conclude by induction that d214l) holds for all integers n ^ 0. 


(215) 

(216) 

(217) 

(218) 

(219) 

( 220 ) 

( 221 ) 

( 222 ) 
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The sequence {ft n ) converges to ft) in L 1 (R) norm and so a subsequence of ft) n converges to ft) pointwise 
almost everywhere. From (12131) and (12141 ) we conclude that 

1^(01 < /3exp(—/x|^|) (223) 

for almost all |£| ^ y/2\. By changing the values of ft) on a set of measure 0 (which does not affect the 
value of R [ft)]), we ensure that (12231) holds for all |£| ^ \/ 2 A. By inserting d213l) into (12231 1 we obtain 

1^(01 < Texp(-//|£|) for all |£| sS V2X, 

which is the conclusion (12051) . 


(224) 


We combine (12231) with Theorem [16] (the application of which is justified since ft < 2 T < 2 A 2 //) to 
conclude that 


IVKOI < exp(-/i|£|) 




A 2 




■ exp 


ft 


47tA 2 /x 


exp 


ft 


47tA 2 


1^1 


(225) 


\ 167T 27T 

for all £ e R. Note the distinction between (1214b and d225b is that the former only holds for £ in the 
interval [—\/2A, v'TA], while the later holds for all £ on the real line. It follows from (1203b that 


1 1 A V 2T 1 

\p 2 'A 

We insert these bounds into (12251) in order to conclude that 


IlHOI < r exp(—a|£|) 




A 


1 


A 7 V 167T 27T 


— exp — exp 


87T 


47tA 


lei 


for all £ e R. We conclude from Lemma |3] that 


exp 


We observe that 

exp 


1 


47rA 

for all £ e R, and that 


lei 


47rA 


I 167T 


Ifl 


A 


sc - 


47T 


exp(l) 


for all £ e 


2 vr 


— exp — exp 


87T 


47tA 


Ifl 




47r 


AT- 


1 


1 


exp 


exp(l) J \l 67 r 27 t 

We combine (1228b . (12291) and d230l) in order to conclude that 


1 

87T 


( 167T 


< 1 . 


2 ^ eXP 


1 

87T 


exp - 


1 


2vrA 


l£l 


fP 

A 


m (. 


i 


i 


A ) \ 167T 27T 


exp 


1 

87T 


exp 


47rA 


1^1 


<: 


ft 2 


for all £ e R. From d210l) and d213b we obtain 


ft 2 ft 2 

A + _ A 

By inserting d231b into d225l) we arrive at 

4T 


which is d206b . 


1^(01 < ( 1 + — ) r exp ( - ( //- 


2 r-r = 2/3-r< (i + yjr. 


l£l for all £ e 


2vrA 


(226) 

(227) 

(228) 

(229) 

(230) 

(231) 

(232) 

(233) 

□ 


Suppose that p is an element of L 1 (R), and that there exist positive real numbers // and T such that 

|p(£)| T exp (-/x |£|) for all £ e R. (234) 

Suppose further that ft) e L 1 (R) is the solution of d!51l) obtained by invoking Theorem [T5l Then the 
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function a b defined by the formula 


°b(x) = 77 - f exp (ixg)ip(g) dt, (235) 

^ J—00 

is a solution of the integral equation (11441) . Moreover, according to Theorem [T7j if A > 2 p " 1 and 
A > 2T, then the Fourier transform of a b (which is, of course, ip) decays faster than any polynomial. In 
this event, cr fo is infinitely differentiable, a b e L 2 (R), and a b e L 2 (R). We record these observations in 
the following theorem. 


Theorem 18. Suppose that there exist real numbers A > 0, T > 0 and p > 0 such that 


A > 2 max 



Suppose also that p is an element of L 1 (R) such that 

|p(f)l < r exp(-p|£|) for all £ e R. 

Then there exists a solution ab e L 2 (R) n C 00 (R) of the integral equation (11441 ) such that 


I<3,(01 < 



T exp (—p|£|) 


for all |£| ^ V2X, 


and 


1^(01 < 



T exp 



for all (eR. 


(236) 

(237) 

(238) 

(239) 


7. Perturbed integral equation 


Suppose that a b is the function obtained by invoking Theorem[[ 8 ]so that 

a b (x) = S [T b [cr b ]] (x) + p(x). (240) 

We rearrange d240b as 

a b (x) = S [T [cr b ]] (x) + p(x) + ( S \T b [a b ]] (x) - S [T [<r fc ]] (x)) (241) 

and define u b via the formula 


u b {x) = S [T b [cj fe ]] (x) — S[T [u b ]] (x) 


(242) 


so that a b is a solution of the perturbed integral equation 


a b (x) = S[T [ab]] (x) + p(x) + v b (x). 


(243) 


From the discussion in Section|4j we see that the phase function arising from a b approximates solutions 
of (jU with accuracy on the order of A~’ ||z^,|j x . However, it is not immediately apparent that T [a b ] is 
defined: the integral 


1 

2A 



y\)<Xb(y ) dy 


(244) 


need not converge absolutely, and without an estimate on the derivative of a b , it is not clear that 

<3,(0 

4A 2 - C 2 ’ 


(245) 


which is formally the Fourier transform of T [uj, defines a tempered distribution (although, in fact, 
it does). It is possible obtain a bound on the derivative of d b by modifying the argument of Section[ 6 l 
That bound can be used show that T [a b ] is defined and to estimate the magnitude of u b . We prefer the 
following, simpler approach to constructing an appropriate perturbation u of the function p. 
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We define the function a via the formula 

?(<£) = VbimO: (246) 

where b(£) is the function used to define the operator 7],. We observe that, unlike T \ a b ], there is 
no difficulty in defining T[a] since a e L 2 (R) and the support of a is contained in (—\/2A, \/2A). 
Moreover, T b \a b \ = T [<r] so that 

(T b (x ) = S [T [cr]] (x) + p(x) for all ieI. (247) 

Rearranging d247b . we obtain 

a(x) = S [T [cr]] (x) + p(x) + v(x) for all ieI, (248) 

where v is defined the formula 

v(x) = a(x) — a b {x ). (249) 


Using (12391) . (1236b and (12491 1 , we conclude that under the hypotheses of Theorem [181 



By combining Theorem [18] with (1250b we arrive at the following theorem. 


(250) 


Theorem 19. Suppose that q e C 00 


x{t) = f Vq(u) 

Jo 

Suppose also that p{x) is defined via the formula 


is strictly positive , and that x(t) is defined by theformida 

t ' _ 

du. 


(251) 


p(x) = 2{t,x}\ 


(252) 


that is, p(x) is twice the Schwarzian derivative of the variable t with respect to the variable x defined via (1251 D . 
Suppose furthermore that there exist positive real numbers A, T and p such that 


A > 2 max 



(253) 


and 


[?(£)[ ^ Texp (—p [£[) for all f e R. (254) 

Then there exist functions v and a in L 2 (R) n C°° (R) such that a is a solution of the nonlinear integral 
equation 

a(x) = S [T[a]](x) + p(x) + u(x), (255) 

p«)is (i + 


2T\ 

—ffj T exp ( /^|Cl) for all |^| < V2A, 


(256) 


and 


a(f ) = 0 for all ]£[ ^ V2X, 


v 




exp (—pX) ■ 


(257) 

(258) 
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The function a b is obtained as the inverse Fourier transform of the limit f of a sequence of fixed point 
iterates {'ip n \ for the equation (11511) . As a consequence of Theorem [2} the functions are elements 
of the space 5(E) of Schwartz functions if p is an element of 5(E). Thus a b is the limit in L l (E) of a 
sequence of Schwartz functions. The following proof of the principal result of this paper, Theoremfl2l 
proceeds by approximating a b using the inverse Fourier transform of an appropriately chosen ip n . 


Proof of Theorem U2l We denote by { p; n } the sequence of fixed point iterates for (11511) generated by the 
function p. That is, ipo is defined by the formula 


^o(0 = p(0 

(259) 

and, for each integer n = 0,1,2,..., V’ n+ \ is defined via the formula 


f>n+ 1(0 = R[i>n] (Oj 

(260) 

where R is as in (|150|). According to Theorem H5l the sequence {'W n } converges in L 1 (R) to a function 
such that 

f’iO = R [V’] (£) for all (eK. 

(261) 

Moreover, as a consequence of Theorem [2] and the assumption that 
element of 5(E). We denote the inverse Fourier transform of d by cr\ r 
a solution of the nonlinear integral equation 

p e 5(E), each of the f> n is an 
According to Theorem Il8l d b is 

d b (x) = S [T b [d b ]] ( x ) + p(x). 

(262) 

We now define a via the formula 


d(£) = HOKO, 

(263) 

where b is given in (l69l). so that 


T b [d b \(0 = T[a] (0- 

(264) 

And we define u via the formula 


*(£) = ^(0 -<rb(€)- 

(265) 

From (|262D, (12641) and (|265|). we conclude that 


d(x) = S[T [<r]] (a;) + p(x) + u(x) 

(266) 

for all x £ E. According to Theorem H9l 


„ T / 4T\ 

l|i>ll ” < v ( 1 + Tj Ktp( -" A) - 

(267) 

For each integer n = 0,1,2,..., we denote by a n the inverse Fourier transform of the function 

Vvi (£)&(£)■ 

(268) 


The function f> n is in 5(E) and b is an element of C'f (R), so the a n are contained in 5(E). We combine 
(12631) and (12681) in order to obtain 


W&n &\\ 




^ J— 00 ^ J — 00 


as n ^ co. 


(269) 
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Since b is supported on (—\/2A, y/2X), 

\\T[a n ]-T[a]\ 


00 ^ 2tt 


< 


I 


47rA 2 


■00 

I 


/—00 
'•OO 


4A 2 - £ 2 


da 


da 


(270) 




4 ^L, IWf)_,S(f)l d( ’ 


from which we conclude that T [cr n ] converges to T \d\ in L x (R). Similarly, if we denote the derivative 
of the function T [f] ( x ) with respect to .x byT[/]' (x), then 


\\T[a n ]'-T [a]'\\ 




i r° 

27T J_ 


-iaHO'tpnio + iaKOHO 


< 


i 


4A 2 - a 2 

r 00 

I |VVi(£) “ 

J — 00 




(271) 


in order to conclude that 

—*• 0 as n —> oo. 


(272) 

(273) 


2\^2tt X J—oo 

from which we conclude that the derivative of T [cr n ] converges to the derivative of T [<r] in L° 

We combine these observations with the definition (I%1) 

|S[T[a n ]]-S[T[a]]|| M 

We rearrange (1266b as 

(T n (x) = S [T [cr n ]] ( X ) + p(x) + n n (x), 

where v n is defined via the formula 

u n (x) = D(x) + (S [T [cr]] (x) - S[T [cr n ]]) (x) + (a n (x) - a(x)) . (274) 

We combine (1274b with (1269b and (1272b in order to conclude that 

K — v Hop —> 0 as n —> oo. (275) 

(276) 


Together (1275b and (1267b imply 


'rail 00 


2p 


4T\ 


< 7T" 1 + T eX P (~^ A ) 


A J 


when n is sufficiently large. Note that the inequality (1267b is strict. 

We have already established that 4> n and a n are elements of 
we observe that 


T[a n m = 


dn(a) 


for all nonnegative integers n. Now 
(£)&(£) (277) 


4A 2 - £ 2 4A 2 - ^ ' 

Since the support of b is contained in (— y/2\, \/2A) and yj n is an element of S'(R), we conclude from 
(1277b that the Fourier transform of T [er n ] — and hence T [er n ] — is an element of S (R) for each non¬ 
negative integer n. Next, we combine this observation with Theorem Q] in order to conclude that 
5 [T \'ip n 11 6 ,5'(R) for all nonnegative integers n. We rearrange (1273b as 

u n {x) = a n (x) — S[T [a n ]] (x) - p(x) (278) 

and observe that all of the functions appearing on the right-hand side of ( 1278b are elements of S(R). 
We conclude that //„ e S(M) for all nonnegative integers n. 

It follows from Formulas (1213b and (1214b . which appear in the proof of Theorem[l8j that 


|^n(S)| < ( 1 + y) r eXP 


(279) 
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for all |^| ^ \/2A and all nonnegative integers n. We conclude from 112731) (12761) and (12791) . and our ob¬ 
servation that a n and v n are elements of 5(R) for all nonnegative integers n that we obtain Theoremll2l 
by letting cr = a n and u = u n for a sufficiently large n. □ 
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